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We investigate the entanglement properties of resonating valence bond states on a two dimensional 
lattice in the presence of dopants that remove electrons from the lattice creating holes. The move- 
ment of the holes generated by the Hubbard Hamiltonian in the regime of strong Coloumb repulsion 
in this setting could be responsible for the phenomenon of high temperature superconductivity as 
hypothesised by Anderson in Science 235, 1196, (1987). We argue that there is a particular density 
of dopants (holes) where the entanglement contained in the lattice attains its maximal value for the 
nearest-neighbour RVB liquid state. 



INTRODUCTION. 

In quantum many-body physics, resonating- valence- 
bond (RVB) states have received a lot of attention due 
to their importance in the description of different phe- 
nomena. They are used to describe the resonance of co- 
valent bonds in organic molecules, behavior of Mott in- 
sulators without long-range antiferromagnetic order [T], 
superconductivity in organic solids [2 , and the recently 
discovered insulator-superconductor transition in boron- 
doped diamond [3]. There are many other applications 
of RVB states (see e.g. [4]). Moreover, RVB states have 
been suggested as a basis for fault-tolerant topological 
quantum computation [5 . 

It was postulated by Anderson in Ref. [Ij that the 
short range nearest-neighbour RVB state (also called 
RVB liquid) might be responsible for the phenomenon of 
high temperature superconductivity. The cuprate super- 
conductors are recognized as doped Mott insulators as in 
the case of Strontium doped Lanthanum cuprate LSCO. 
The copper oxide planes are described by the one-band 
Hubbard model with a strong, on-site repulsion U. The 
RVB state on a square lattice for the Hubbard Hamil- 
tonian with strong Coloumb repulsion was proposed by 
Anderson to be the Mott insulator phase of the system, 
i.e., the pure Lanthanum cuprate is in an RVB state. By 
introducing dopants one removes some electrons from the 
lattice creating "holes", i.e., unoccupied sites on the lat- 
tice. As pointed out by Anderson, it would take a finite 
concentration of dopants to metallize the material and 
make it a superconductor as initially the dopants will be 
screened by the bound quasiparticles. The holes start 
to hop from one site to another and their motion is un- 
inhibited leading to their delocalization over the whole 
lattice, which can be interpreted as a persistent current. 
The movement of the holes is resistant to thermal noise 
and depends on the density /i of dopants (holes) reach- 
ing its maximum somewhere between 10% to 15% (the 
so-called 1/8 anomaly) of holes on the lattice. It is no- 
table that there does not seem to be a clear consensus on 
the origin of this anomaly. 



Indeed, the ground state of the 2D Hubbard model 
with doping is still unknown. Numerical simulations in- 
dicate that the RVB scenario is the right one for cou- 
pled plaquettes and ladders [6 , and recently experiments 
have been proposed to test the RVB scenario in fermionic 
atoms in 2D optical lattices [7J. These experiments pro- 
pose methods to increase the inter-ladder coupling to 
check if the RVB state on ladders is adiabatically con- 
nected to the Hubbard model ground state on the square 
lattice, which would provide an experimental test of the 
RVB theory. In this paper, we investigate the entangle- 
ment properties of the RVB states in small-sized laddders 
and plaquettes and speculate on the behavior in the ther- 
modynamic limit. 

The entanglement properties of RVB states without 
dopants on many dimensional lattices has been investi- 
gated in [8 . The main conclusion of the Ref. [8 is that 
such states can only have (if any) a very small amount 
of bipartite entanglement between any two sites on the 
lattice but they are always genuinely multi-partite entan- 
gled. 

In view of this, it is interesting to see how multipartite 
entanglement of RVB states depends on the density of 
holes. In particular, we are interested if changes in the 
amount of entanglement correspond to the experimental 
observation of the maximal Tc superconductivity in the 
hole density window 0.1 < ji < 0.15. Intuitively, one 
would expect entanglement to be larger in this window as 
well. This is because entanglement and supercurrent are 
both related to the existence of correlations and therefore 
their peaks should be related i.e. a larger entanglement 
would be more robust to increase in temperature. 

To address these questions we first define the RVB 
states on a 2L site lattice with 2n holes and an appropri- 
ate entanglement measure. Subsequently, we investigate 
analytically lattices up to 24 sites with open boundary 
conditions. Based on the obtained results we conjecture 
that in the thermodynamic limit of L ^ oo the amount 
of entanglement reaches its maximum for some critical 
density of holes ficn where < jicr < 1- 
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FORMULATION OF THE PROBLEM 

Let us consider a two dimensional lattice with open 
boundary conditions consisting of 2L sites, which is a 
union of two sub-lattices A and B in such a way that 
any site belonging to the sub-lattice A{B) has all its 
nearest neighbours belonging to the sub-lattice B{A). 
We define a dimer between sites a G A and b G B 
as a singlet state \Sab) = ^(|0)a|l)6 - |l)a|0)6). A 
dimer covering is defined as a tensor product of dimers 
I^[(ai5i),(a262),...,(az.6z.)]) = ^^Lil^a^bJ. Here the set 
of pairs [(ai^i), (a262), . . . , (a^^L)] represents a partic- 
ular way of joining neighbouring sites of the two sub- 
lattices with singlets. The number of such sets for the 
case of square lattices with open boundary conditions 

is, from [2], m, given by [1^(4 cos2(2^) + 

4cos^(2^^)). For periodic boundary conditions in the 
square lattice this number is known to be (exp (^))^ ^ 
(1.791)^, where G = T.n^oi-'^T / i^^ + 1)^ is Catalan's 
constant. 

The RVB liquid is defined as 

^ 77^ I^[(ai&i),(a262),...,(az.6i.)]), (1) 

{aibi)...{aLbL) 

where R is the normalization constant and the summa- 
tion extends over all possible dimer coverings. R can 
be in principle calculated using the techniques of the so- 
called 'random loop soup' [TT]. In practice, the problem 
of counting is analytically intractable and even with the 
help of numerics one cannot compute R for large L. 

Let us now denote an equal superposition of dimer cov- 
erings of the lattice where two arbitrary sites and bj 
are unoccupied by |A(a.5^.)). Please note that we do not 
allow two sites belonging to the same sub-lattice to be 
unoccupied as this will preclude the possibility of cover- 
ing the rest of the lattice with nearest-neighbour dimers. 
This guarantees that the occupied part of the lattice is in 
the RVB liquid state with the number of dimer coverings 
depending on the positions of the holes. 

In general, if there are 2n holes one can define 
in the same manner the state \^[{aibi),ia2b2),...MM]) ^ 
i.e., the superposition of all possible coverings of 
the part of the lattice excluding the empty sites 
(ai^i), (a262), . . . , {dnbn)' Note that if n is too large it 
may not be possible to cover the part of the lattice with- 
out holes by dimers. A simple example is when 2n = L. 
In this situation if every second site is empty (chess board 
configuration), the remaining part of the lattice cannot 
be covered by dimers. We will address this issue further 
on. 

The final steady state of the system, i.e., the state 
after the system of the lattice and dopants has reached 



an equilibrium can be written then as 

|A2„) = Yl 

iai-^bj^),...{a^^bjJ 

where the probability distribution of the holes 
p{ai^ ^bj^^ . . . , bj^) depends on their detailed dy- 
namics. 

We would like to be very clear about the meaning of 
I A2n), namely that, it may not be an accurate description 
of the state of the lattice of a real doped superconductor. 
First of all, the RVB theory is one of many theories of 
the discussed phenomenon [12 . Secondly, even within 
the RVB theory itself it is not clear how the state with 2n 
holes looks like [13^. For instance, one cannot exclude the 
possibility that long-range dimers will appear in |A2n)- 
However, it seems reasonable to assume that for a small 
amount of holes |A2n) is an acceptable choice. 

It is not clear how to quantify multiparty entangle- 
ment in the |A2n) state. The problem arises from the 
fact that a commonly accepted definition of multiparty 
entanglement is that any bipartition of the considered 
many particle quantum state must be entangled. How- 
ever, in our case any bipartition of the lattice (equiva- 
lently the bipartition of the state |A2n)) leads to a state 
with a variable number of particles on each site of the bi- 
partition. According to super-selection rules one cannot 
observe a superposition of states with different number 
of particles which considerably complicates the task 
of quantifying entanglement in RVB states with holes. 

As a measure of the amount of non-classical correla- 
tions in the lattice we take the geometric measure of en- 
tanglement [15 , which is generalized to the multi-partite 
case in a straight-forward manner [16]. For pure states, 
the measure is given by the ^ -based logarithm of the 
squared modulo of the overlap between the state, and 
the separable state closest to it 

i?(|^/;)) = -max|log2|(V'|V«ep)|', (3) 

where {tpsep) is a separable state. In the case of our lat- 
tice, however, we deal with the subtle matter of super- 
selection rules, because the closest product state to the 
RVB state could involve forbidden local superpositions 
of a hole and an electron. For this reason, we use the 
average geometric measure. The averaging is done over 
all possible locations of the holes. 

More precisely, we define an average geometric measure 
of entanglement on the state | A2n) as 

E{2n) = 

{ai^bj^),...{ai^bj^) 

p{ai, ,bj,,... ai^,bjjE{\A^^a^b^),...,iaM]))^ (4) 
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where 

^(|A[(ai6i),(a262),...,(a^6^)])) = 

-21og2 max \{^sep\^[{a^h^),{a2h2)....MM])\' (^) 

The maximum is taken over ah fuhy separable states on 
the part of the lattice without holes and the physical 
meaning of E(2n) is clear; it is the average amount of 
entanglement one gets after locating the position of the 
holes on the lattice. 



RVB state without holes: ^(0). 

First we consider the RVB state without holes. We 
already know from the Ref. [8 that it contains negligible 
two-site entanglement but it is genuinely multi-party en- 
tangled. Here we calculate ^(0) that will serve us later 
as a basis for comparison with E(2n). Additionally, the 
same method of calculation will be used for n > 0. 

We have 

^(0) = -2 max log^ |(^.ep|A)| = -21og2 ^, (6) 

where C is the number of dimer coverings for |A), i.e., 

C = U^iU&'icosH^) + 4cos2(5^)) for open 
boundary conditions and C = (exp(^))^ for periodic 
boundary conditions in the square lattice. This result 
can be argued as follows. The state |A) always has 
an " anti- ferromagnetic" term of the form |0101...01) 
or 1 1010 ... 10) with its coefficient equal to Natu- 
rally every other term has smaller coefficient, because in 
the superposition of all coverings only the antiferromag- 
netic terms add up. Thus, there is a fully separable state 
IV^sep) = 1 0101 ... 01) for which the modulus of the scalar 
product with |A) equals ^. However, the only fully sep- 
arable state with equal number of zeros and ones is of 

the form \x1X2 . . . Xn) with Xi + X2 H \-Xn = which 

means that = iV^sI?^'''^^). 




FIG. 1: E{2n) for Ladder lattices 3x2, 4x2, 10x2. Note 
the appearance of the peak at 2 holes for lattices larger than 
6x2 

|A[(ai6i),(a262),...,(a.6.)])- Therefore, we have 

E(2n)= ^ p(ai,6i,...an,6n) 

^Qg^ f ^[(aibi),(a262),-,(a.6^)] \ 
V^[(ai6i),(a262),...,(a^6^)]/ 

The main difficulty in the above formula is that it is an 
A^P-complete problem to calculate C^^alhl),{a2h2),...,{ar^hr^)] 
and ^[(al6l),(a262),•••,(ar^6r^)]• ^^^^^ stcms from the known 
result in theoretical computer science [17 that count- 
ing dimer coverings of a planar lattice is a polynomial- 
time computable problem, whereas counting monomer- 
dimer arrangements on the two-dimensional lattice is an 
A^P-complete problem. In our language, the holes cor- 
respond to the monomers and the singlet pairs corre- 
spond to the dimers. Thus, finding C'[(ai6i),(a262),...,(a^6^)] 
is an A^P-complete problem. Moreover, calculating 
^[(ai6i),(a262),...,(a^6^)] IS at bcst polyuomial-timc re- 
ducible to finding C[^arhr),{a2h2)^...MM]' 



RVB with 2n holes: E{2n). 



RESULTS 



To compute entanglement in this case it suffices to find 
the maximal overlap between | A [(^^ 5^)^(^262 ),•••, (a. &r.)]) 
and a fully separable state for every possible set 
of pairs (ai^i), (a262), . . . , (a^^n)- From the previ- 
ous considerations we know that the maximum is 
reached for an anti-ferromagnetic separable state and it 

^g^^g ^[(aib,),(a2.2),-,(a.b.)] ^Yi^^^ CuaM (a,b2) (a b )] 

^[(aibi),(a2b2),...,(a^b^)] ' Ll^l J ^^2 02 j , • • • , l^r^ 0^ j J 

is the number of the coverings of the initial lattice 
with the sites (ai, 61), (a2, 62), • • • , {dn, ^n) removed and 
^[(al6l),(a262),•••,(or^6r^)] thc uormalizatiou of the state 



In this section we present the main results of the paper. 

We have analytically computed entanglement E{2n) 
using the method described above for ladders and pla- 
quettes up to size 6x4 (ladder lattices of size 3 x 2, 4 x 
2, . . . , 10 X 2 and rectangular lattices 4x3,4x4, 5x4 
and 6x4). The results, shown in the Fig. [1] for lad- 
ders and Fig [2] for rectangular lattices, clearly show that 
E{2n) reaches the maximum at a certain hole density as 
the size of the lattice increases. Moreover, the maximum 
occurs at a low concentration of holes. 
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FIG. 2: E{2n) for Rectangular lattices 6x4, 5x4, 4x4, 4x3. 
Note the appearance of the peak at 4 holes in 5x4 and 6x4. 
The position of the peak has shifted to the right (higher den- 
sity) compared to that in Fig [1]. 



We now elucidate certain significant points in the cal- 
culations leading to the graphs in Figs [1] and [2]. As an 
illustrative example, let us consider the lattice of size 6x4 
with four holes, at which point the peak occurs in this 
structure. It is clear that the 4 holes can be in one of 
X C^^ positions, where denotes the binomial coef- 
ficient. Hence, one has to average over the entanglement 
found in each of these cases to find ^(4) for this lattice. 
However, in this calculation we omit the pathological po- 
sitions of the holes in which a single site is surrounded on 
all sides by holes as in Fig. [3], in which case the rest of 
the lattice is unable to form a short-range RVB structure. 
Long-range dimers between sites belonging to the same 
sublattice would be needed to fill the lattice in such a sit- 
uation and we omit the corrections accruing due to these. 
In any case, neglecting these situations cannot substan- 
tially alter the behaviour of entanglement because they 
occur with the probability of 4(C^''^/^^"^ - I) / {C^^^^'^^f 
for four holes in an a x 6 lattice. Note that the probability 
of occurrence of these situations is zero for two holes so 
that the positive gradient at the beginning of the curve 
is maintained and the existence of the peak is assured. 
Thus, we conclude that the peak in the graph is main- 
tained even when these situations involving long-range 
dimers are taken into consideration. 

It is seen from Fig [2] that the initial gradient of the 
curve increases with the size of the lattice. If the trend 
continues for larger lattices, one might expect the peak 
to shift to the right and converge to a particular concen- 
tration in the thermodynamic limit. 
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FIG. 3: Typical pathological case which is omitted from short- 
range RVB calculations. These cases do not affect the behav- 
ior of entanglement as explained in the text. 

CONCLUSIONS 

Our data strongly indicates that the average multi- 
party entanglement quantified by E(2n) reaches a max- 
imum for some critical density of holes in the thermo- 
dynamical limit. Although we are unable to predict the 
exact value of this critical density due to the hard com- 
putational nature of the problem, we conjecture that it 
is located in the region of low density of holes. 

A possible way to get some additional information 
about the location of the maximum in the thermody- 
namic limit would be to translate the problem in graph 
theoretic language in the following way [l^ • It is appar- 
ent that each dimer covering of the RVB structure can be 
represented as a balanced bipartite graph with the sub- 
lattices A and B forming the two vertex sets, there being 
L edges in a 2L site lattice. In such a graph, there is no 
path between any two sites belonging to the same sub- 
lattice, and the matching number of the graph is equal 
to L. These conditions can be extended to the case of the 
lattice with holes as well. 

The problem of finding C then translates to the equiv- 
alent problem of finding the number of graphs with ad- 
jacency matrices having the following property. Each 
adjacency matrix is a sparse matrix of size 2L by 2L, the 
mth row has a 1 at one of upto four possible positions, 
these corresponding to the sites adjacent to site m in the 
lattice; and each row and each column has only one non- 
zero entry. The number of such matrices then equals the 
number of coverings. Algorithms for approximating C 
have been devised as enumerating the number of perfect 
matchings in such graphs in [19 . 

The problem of finding R can also be broken down into 
the equivalent problem of finding the number of degener- 
ate and non-degenerate loops in every distinct superposi- 
tion of two coverings, as in [20]. To do this, we add all the 
adjacency matrices for that lattice, two at a time, keep- 
ing only the distinct results. The number of degenerate 
loops (dl) in each such superposition is then simply half 
the number of twos in the resulting matrix. The number 
of nondegenerate loops (ndl) is equal to the total num- 
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FIG. 4: Graph theoretic formulation: 1(A), 2 (A)- Sample cov- 
erings 1 and 2 of 4 X 4 lattice with two holes; 1(B),2(B)- 
Corresponding bipartite graphs; 1(C), 2(C) Corresponding ad- 
jacency matrices to the two coverings 




(A) Superposition of Coverings 1 and 2 

(B) Graph 

(C) Adjacency Matrix of Resulting Graph; 
Number of degenerate loops = 2 ; 
Number of non-degenerate loops =1 



□ 2000000000000 
20000000000000 
00020000000000 
00200000000000 
OOOOOIOOOOOIQQ 
OOOOIOIOOOOOQQ 
00000101000000 
00000010100000 
00000001010000 
00000000101000 
00000000010001 
00001000000010 
00000000000101 
00000000001010 



FIG. 5: Graph theoretic formulation: Coverings 1 and 2 su- 
perimposed on each other 



ber of cycles in the graph for which the resulting matrix 
forms the adjacency matrix, minus the number of degen- 
erate loops. R is then found from the neat formula, [20] 



R 



E 

superpositions 



2dl ^ 4r 



(8) 



However, there do not seem to be good approximation 
schemes to bind R using the above method. 

Figs [4] and [5] illustrate the method in graphic de- 
tail showing two coverings for a 4 x 4 lattice with two 
holes. In Fig [4] 1(A), the red circles indicate sublat- 
tice A sites and blue circles indicate those belonging to 
sublattice B. Shaded circles indicate holes in the lattice 
and singlets are represented as lines between sites. The 



bipartite graph corresponding to the state is shown to 
its right in 1 (B). The graph has bipartition (A,B) with 
A = (1,3,5,7,9,11,13) and B = (2,4,6,8,10,12,14). 
Since, \A\ = \B\ = L, the graph is balanced. Edges of 
the graph connect vertex set A to B such that an edge 
connects a vertex to only one of its nearest neighbouring 
vertices on the lattice. Since there are L such edges, the 
size of the maximum matching of the graph is equal to L. 
The adjacency matrix for this graph is shown alongside 
in 1(C). To find R, we would need to superimpose the two 
coverings on each other as shown in Fig [5]. The result- 
ing graph and its adjacency matrix (the sum of the two 
adjacency matrices in Fig [4]) are shown alongside. The 
number of degenerate and non-degenerate loops can then 
be calculated from the number of cycles in the graph. 

It is hoped that with these methods and by experimen- 
tal observations as suggested in [7] the entanglement vs 
hole density curve can be constructed in the thermody- 
namic limit. This might throw more light on the question 
of whether multipartite entanglement, defined in this av- 
erage geometric sense could indicate the occurrence of 
the quantum phase transition. 
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